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Abstract A profound change occurs in the stability of quantized vortices in exter- 
nally applied flow of superfiuid 3 He-B at temperatures < 0.6 T c , owing to the rapidly 
decreasing damping in vortex motion with decreasing temperature. At low damping 
an evolving vortex may become unstable and generate a new independent vortex loop. 
This single-vortex instability is the generic precursor of turbulence. We investigate the 
instability with non-invasive NMR measurements on a rotating cylindrical sample in 
the intermediate temperature regime (0.3 - 0.6) T c . From comparisons with numerical 
calculations we interpret that the instability occurs at the container wall, when the 
vortex end moves along the wall in applied flow. 

Keywords quantized vortex • vortex formation ■ vortex instability • vortex dynamics ■ 
turbulence • transition to turbulence • mutual friction • rotating flow • vortex remanence 

PACS 67.57.Fg ■ 47.32.-y • 67.40.Vs 
1 INTRODUCTION 

Turbulence in a rotating fluid is a most frequently encountered phenomenon, ranging 
from astrophysical and planetary scales to engineering problems. Turbulence in a ro- 
tating superfiuid has not been extensively investigated, although it is a simpler form 
of turbulence from that in viscous fluids and might provide a shortcut to new under- 
standing. At temperatures below 2.5 mK, uniform rotation is one of the few technically 
feasible means of generating flow on the global scale. The main reason for studying 
superfiuid 3 He-B at these temperatures is to learn about the influence of a strongly 
temperature dependent vortex damping [1 , or mutual friction dissipation a(T). With 
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decreasing temperature a sudden onset of turbulence is seen in applied flow. It is this 
onset and the mechanisms behind it that are described in this report. An instability 
of a single vortex evolving in applied flow is at the root of the onset. We study the 
instability using a new measuring technique, the injection of a seed vortex in rotat- 
ing vortex-free counterflow, and monitor the subsequent evolution of vorticity with 
noninvasive NMR measurement |2]. 

Spin up of the superfluid component: — In uniform rotation at an angular 
velocity fi the induced applied flow is the superfluid counterflow (cf) of the normal and 
superfluid fractions, the relative velocity v = v n — v s . The viscous normal fraction (with 
olive oil-like viscosity) is practically always in solid-body rotation with the container, 
v n = n x r, while the velocity v s of the superfluid fraction is produced by the combined 
flow field from all quantized vortex lines and from the flow caused by the boundary 
conditions on the container walls. Superflow is characterized by a finite critical velocity 
Vc,exp, at which vortex formation starts. If rotation is started at temperatures below 
T c , then the vortex-free Landau state is formed first, where v s = (in the stationary 
laboratory frame). When the cf velocity reaches the critical value u c ,cxp (somewhere 
at a rough spot on the container wall), a vortex is formed, which then evolves from 
a micron-size loop to a rectilinear line vortex in the center of the cylinder, aligned 
parallel to the rotation axis. Thereby the maximum cf velocity drops below the critical 
limit, v(R) < i> c ,cxp- If the external rotation drive fl is continuously increased, then 
the process repeats and a central cluster of rectilinear line vortices is formed which is 
contained and isolated by an annular layer of vortex-free cf from the cylinder wall. The 
maximum cf velocity in the vortex-free flow is at the cylinder wall, with v(R) ~ u c ,cxp- 
Within the vortex cluster the global cf vanishes (on length scales exceeding the inter- 
vortex distance d v ) and v s is on an average in co-rotation with v n . 

This is the linear well-behaved vortex formation process above 0.6 T c at high mu- 
tual friction damping, where it provides the "spin up" of the superfluid component to 
co-flow with the normal fraction. In macroscopic flow geometries such "quasi-intrinsic" 
vortex formation can be observed in superfluid 3 He-B [3J, where (depending on surface 
roughness of the cylinder wall) the experimental « c ,oxp is perhaps only an order of 
magnitude smaller than the theoretical v c bulk; while in 4 He-II this has been demon- 
strated only with flow through orifices of < 1 ^tm size. Typically in a smooth-walled 
cylinder u c ,exp ~ 1 cm/s at about 0.7 T c from where it decreases approximately as 
^c,exp OC - T/T c towards T c . Although this is an order of magnitude smaller than 
the bulk liquid critical velocity w Ci bulk ~ 11 cm/s (at P = 29 bar liquid pressure), 
vortex-free cf can thus be maintained in metastable state up to v(R) < u c ,cxp, which 
substantially influences the NMR absorption response. 

In 3 He-B mutual friction decreases towards low temperatures almost exponentially. 
Thus the transition from the laminar to the turbulent flow regime at 0.6 T c is abrupt and 
sudden pQ. Below 0.6 T c mutual friction dissipation becomes sufficiently small so that 
an evolving vortex may become unstable and generates a new vortex loop which then 
starts expanding independently. This can happen repeatedly at constant rotation [I]. 
The experimental signature of the instability is a sudden burst of turbulence, when the 
density of evolving vortices grows sufficiently so that they start interacting turbulently 
in the bulk volume. To characterize this chain of events, we need to know how the 
first new loop is generated from a vortex which is evolving in applied flow. This is 
investigated here using the scheme outlined in Fig. [T] 

The mechanism by which the superfluid fraction is set into rotation is often called 
"spin up" of the superfluid component [5] . A sequence of processes is involved in spin 
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Fig. 1 (Color online) Principle of measurements on dynamic vortex generation. Well separated 
isolated seed vortices are introduced in rotating vortex-free cf. The initial high-energy state of 
vortex-free flow may then relax to the equilibrium vortex state via vortex generation processes 
which become possible at temperatures below the hydrodynamic transition at 0.59 T c (at 
29 bar liquid pressure). The instability of a single evolving seed vortex is the first step. It is 
followed by a turbulent burst which is started when the density of evolving vortices is sufficient. 
The probability of the combined process depends on the dynamic mutual friction parameter 
£ = (1 — a') J a which is shown on the top together with its range of values in superfluid 4 He. 
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up, one of them being the formation of new vortices. Ideally, our approach here is to 
study spin-up at constant rotation ft, by injecting a seed vortex in vortex-free flow. 
On the macroscopic scale, superfiuids often mimic the behavior of viscous liquids [H]. 
However, as we shall see below, our superfluid spin up is rather different from the 
viscous flow patterns which evolve in the spin up of classical liquids .7 . Owing to its 
high viscosity, the normal component is during the spin up in a state of laminar flow. 
This is an important simplification over the complicated coupled turbulent spin up of 
the normal and superfluid fractions which is observed in 4 He-II at higher temperatures 
i. 

Generic properties of single-vortex instability: — The single- vortex instability 
below 0.6 T c and the phenomena which it starts in a rotating column with circular cross 
section have been described in Ref. [8]. Here we repeat some central features in the 
light of Figs.[2]and|3} before we turn to a closer characterization of the instability itself. 
One type of measurement on evolving vortices in a rotating cylinder is illustrated in 
Fig. [2] This numerical calculation with the vortex filament method [9] (see Sec. [3| 
monitors remanent vortices [10] in a rotating cylinder (with radius R and length L at 
constant J?), in a situation where the vortex instability does not occur. The purpose 
is to illustrate the motion of remanent vortices while they evolve from short curved 
vortices to rectilinear lines. The practical outcome from such a measurement is analyzed 
in Fig. [3] as a function of temperature. 

A measurement of the kind depicted in Fig. [2] proceeds as follows: As shown in the 
inset, an equilibrium vortex state is first decelerated to zero rotation. After a waiting 
time At at zero rotation, some vortices have not yet managed to annihilate. Their 
annihilation time is governed by mutual friction damping a(T) and increases rapidly 
with decreasing temperature. When rotation is next suddenly increased from zero to a 
steady value J?f, the remaining vortices start expanding towards their stable state as 
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Fig. 2 (Color online) Numerical calculation of vortex motions in a rotating cylinder: (t < 0) 
Initial state with 22 vortices at O.lrad/s rotation. The vortices have been artificially tilted 
by displacing their end points uniformly by 1 mm at both end plates of the cylinder, to break 
cylindrical symmetry. Rotation is then abruptly reduced to zero, to allow vortices to annihilate. 
(t < 600 s) After a waiting period At = 600 s, 12 remanent vortices remain which are here 
shown at Q = 0. Rotation is then increased to Qf = 0.5rad/s (t > 600 s) and the 12 remnants 
start evolving towards rectilinear lines. This requires that the vortex ends on the cylindrical 
wall travel in spiral motion to the respective end plates. The mutual friction parameters are 
a = 0.18 and a' = 0.16 |13| . In the picture the radial lengths have been expanded by two, 
compared to axial distances. 



rectilinear lines. This motion proceeds such that the vortex ends on the cylindrical wall 
travel along a spiral trajectory. During such helical motion on outer evolving vortex is 
wound around the straighter vortices in the center |11| . The spiral motion is evident 
on the far right where we see a cluster of helically twisted vortices in the center of the 
cylinder. This state is still evolving, since ultimately also the helical twist relaxes to 
rectilinear lines, when the vortex ends slowly slide along the end plates of the cylindrical 
container [12] . 

Using the measuring routine in Fig. [2j we can find the onset temperature of the 
single vortex instability. This is done by measuring the number of rectilinear vortex lines 
in the final state at different temperatures. A rough value for the onset temperature 
can be located easily: (1) Well above T on the number of vortices does not change during 
the evolution of the remnants and thus in the final state there are only few vortices. (2) 
Well below T on , in contrast, the instability occurs always and is followed by a turbulent 
burst which starts the evolution towards the equilibrium vortex state [8]. These two 
types of final states correspond to very different NMR absorption spectra. In Fig.[3]the 
distribution of these two final states is plotted in the immediate vicinity of Ton- Each 
data point corresponds to an independent measurement of the final state, after the 
remanent vortices have expanded in the applied cf at constant external conditions (f2f , 
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Fig. 3 (Color online) Measurements on the onset temperature T on of turbulence. The mea- 
surement starts from an initial state which is obtained by decelerating an equilibrium vortex 
state at 1.7rad/s to Q = at a rate 0.01 rad/s 2 . The remaining vortices are left to annihilate 
|10| for a period At at Q = 0. Rotation is then increased to i?f at a rate 0.02 rad/s 2 . When 
all transients have decayed the number of vortices is measured in the final steady state at Qf. 
The result is plotted as a function of temperature with 30 - 40 data points per panel. The 
solid curve is a gaussian fit which represents the probability for turbulence with a half width 
ctt = 0.02 T c and centered around T on . Comparing results in the two panels for At = 20 min 
and 2 min, we see that T on decreases with increasing At, since the number, average size, and 
density of remnants is reduced as At increases. Parameters: R = 3 mm, L = 110 mm, and 
P = 29.0 bar. 

T, and P) and all transients have relaxed. Only in the temperature regime T ~ T on the 
final state is unpredictable. By fitting the measurements to the normal distribution we 
get the probability of starting the turbulent burst, the quantity plotted on the vertical 
scale in Fig. [3] 

Two cases are compared in Fig. [3] on the right the waiting time at zero rotation is 
At = 2 min, while on the left it is At = 20 min. This turns out to yield different values 
for Ton- With At — 2 min, the number of remanent vortices A/j is about ten times 
larger (approximately 60 at about 0.53 T c ) than with At — 20 min (approximately 
10 at about 0.47 T c ) [10| . Because of the lower seed vortex density in the latter case, 
the onset moves from 0.53 T" c to 0.47T C . Similarly, if Qf is reduced from its value of 
0.7rad/s in Fig. |3j the onset moves to a lower temperature. This means that, to start 
turbulence, a sufficiently low mutual friction dissipation a(T) is the most important 
condition, but what also matters is the velocity of the applied cf and the initial number, 
configuration, and density of seed vortices. 

The measurements in Fig. [3] are not conducted with single- vortex resolution: The 
number of rectilinear vortices N in the final state is determined within ±10 lines. While 
the number of seed vortices is small, in the equilibrium vortex state it is N C q ~ 640, 
and thus these two situations are easily distinguished. In the onset regime T ~ T on , 
both failed and successful attempts for a transition to turbulence occur. In the failed 
attempts no increase in the number of vortices is measured, which means that only 
few new vortices are needed before turbulence manages to switch on. Secondly, all 
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final states in Fig. [3] are either equilibrium vortex states or states with essentially no 
new vortices. This suggests that when turbulence is switched on, a surplus of vortices is 
created in the turbulent burst. Simultaneously the polarization of the vortices along the 
rotation axis grows to > 90 % and the number of vortices adjusts itself approximately 
to that of the equilibrium vortex state. 

It needs to be pointed out that the generation of new vortices in Fig. [3] is very 
different from regular vortex formation at temperatures above 0.6 T c . The vortices in 
Fig.[3]are created at constant rotation i?f at a low cf flow velocity < 2mm/s relatively 
evenly along the entire column [3]. If turbulence is started, then the macroscopic cf 
velocity ultimately drops close to zero. This means that vortex generation proceeds 
until completion (with the exception of a few rare examples which have been observed 
only in the onset regime, T w Ton)' In contrast, in regular vortex formation above 
0.6 T c at a rough spot on the cylinder wall (or at a wall defect called a "vortex mill" 
|14p. the flow becomes sub-critical after the first new vortex. A second vortex is not 
created (in the ideal situation [3]), unless Q is again increased. 

Thus to summarize, in contrast to the quasi-intrinsic spin-up process above 0.6 T c , 
the single- vortex instability as the precursor to turbulence has unmistakable features, 
which become evident when one cools down into the intermediate temperature regime, 
from 0.6 to 0.3 T c . Here the quasiparticle mean free path £ < 50/mi is still smaller 
than the typical inter- vortex distance d v ~ 0.2 mm and much less than the sample 
size R = 3 mm. In the following we focus on the single- vortex instability and the 
immediately following turbulent burst. 

2 MEASUREMENTS ON SINGLE- VORTEX INSTABILITY 

Rotating flow states: — Controlled seed vortex injection and the calibration of mea- 
sured NMR signals requires stability and reproducibility of different rotating flow 
states. Above 0.6 T c one can experimentally prepare a state with any number of recti- 
linear vortex lines in the central cluster up to the equilibrium number: N < JV q. These 
are called (i) the vortex-free state [N = 0), (ii) a metastable vortex cluster (JV < JV e q), 
and (iii) the equilibrium vortex state N ss A'eq. When a vortex is formed as a small 
loop, it expands in spiral motion to a rectilinear line and becomes part of the central 
cluster, as seen in Fig. [2] Thereby the radius R of the cluster increases and the cf ve- 
locity is reduced: v — Or — nN/(2irr). Here k = h/(2m^) is the superfluid circulation 
quantum, N is the number of vortices in the cluster, and R < r < R. The cluster 
reaches its maximum radius in the equilibrium vortex state, where R = R — d oq and 
deq > d v ~ \J k/ (2,1?) is the equilibrium width of the vortex-free annulus between the 
cluster and the cylinder wall. Independently of the number of vortices N, the maximum 
applied flow is at the cylindrical boundary. At higher J? (when there is no annihilation 
barrier [T3]) the equilibrium vortex state is also the state with the maximum number 
of vortices (at a given value of Q in stable conditions). 

The precondition for generating these rotating states is a sufficiently high and stable 
critical velocity of vortex formation f c ,exp [3|. Since surface roughness reduces i> c ,cxp, 
smooth and clean container walls are important. Thus v c ,exp is container dependent 
and varies even from one cool down to the next, presumably owing to frozen residual 
gas crystallites on the walls. Since the vortex number is reliably conserved during vortex 
formation processes only above 0.6 T c , a particular rotating state generally has to be 
formed at high temperatures, but can then be cooled to low temperatures at constant 



Fig. 4 (Color online) Sketch of equilibrium vortex state in a slightly tilted cylinder with two 
vortices connecting to the cylindrical side wall. 

rotation. In particular, to have a rotating vortex-free sample at low temperatures, the 
procedure is to warm up first above 0.6 T c , where the annihilation of remnants is rapid 
|10j . The sample is kept there at zero rotation for 10 to 20min and is then cooled down 
in rotation (with Q < v c ,cxp/R). 

The equilibrium vortex state is a particular case. Ideally, when the sample and 
rotation axes are perfectly aligned, all vortices in the equilibrium cluster are rectilinear. 
In practice, some misalignment is unavoidable, typically of order ~ 1°. The result of 
this is that some of the outermost vortices are then curved and attached at one end 
to the cylindrical wall (Fig. |4j). When J? is increased, these vortices start to spiral 
towards rectilinear lines. Below 0.6 T c this evolution may become unstable and leads 
to a turbulent burst, as discussed in the context of Fig. [3] Thus the equilibrium vortex 
state can be used for measurements on the single-vortex instability in a similar manner 
as remanent vortices. 

Seed vortex injection: — By introducing seed vortices in applied flow by exter- 
nally controlled means, we can monitor vortex evolution as a function of time at con- 
stant rotation. In rotating flow any vortex which is not rectilinear can be used as seed 
vortex, i.e. a small vortex loop or any section of a longer vortex whose configuration 
changes appreciably as a function of time. The active section of the seed vortex is 
outside the vortex cluster in the counterflow region, where the macroscopic cf velocity 
v(r) 7^ 0, and is curved, with one or both ends connected to the cylindrical side wall. A 
number of different methods exist to start seed vortex evolution [161 [2]. One example 
is the use of remanent vortices (Figs. [2] and[3]l and a second the curved peripheral 
vortices of the equilibrium vortex state (Fig. 

An efficient injection method is based on the superfluid Kelvin-Helmholtz insta- 
bility 17 of the phase boundary between the A and B phases of superfluid 3 He. In 
this instability a tightly packed bundle of several approximately parallel vortex loops 
escapes across the AB interface into the vortex- free B-phase flow |18| . This allows 
immediate inter- vortex interactions and starts the turbulent burst. There are several 
further important points to note. First, the burst follows instantaneously the KH in- 
stability and no precursory vortex generation via the single vortex instability is ob- 
served before the burst. Secondly, the burst can be localized to the immediate vicinity 
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of the AB phase boundary [2]. Thirdly, the KH transition to turbulence is indepen- 
dent of the applied flow velocity v(f2, N, R) [I], unlike other turbulent bursts which 
require the single- vortex instability as precursor [19] . Most importantly however, com- 
pared to other injection methods, KH injection gives the highest onset temperature 
for turbulence. This means that the turbulent burst can happen at a higher value of 
vortex damping a(T) than where the single-vortex instability becomes possible. In KH 
injection the transition temperature to turbulence displays a typical narrow normal 
distribution [19] which is similar to that in Fig. [3] (except for the value of T on which 
is higher). This similarity suggests that plots of the transition to turbulence, like that 
in Fig. [3] describe the transition probability in a situation when enough vortices have 
already been created by the precursor mechanism so that turbulence can switch on. 
For these reasons KH injection is at present the process which is believed to identify 
most clearly the hydrodynamic phase transition between regular and turbulent vortex 
dynamics. 

Finally, we note that the capture reaction of a thermal neutron in 3 He-B releases 
vortex rings in a predictable manner into vortex- free cf [50]. Vortex formation is here 
explained to happen via the Kibble- Zurek mechanism [21], during a sudden tempera- 
ture quench from the normal phase into the superfluid. This process can be reliably 
adjusted to emit only one single vortex ring into vortex-free B-phase flow. The telltale 
observation is that even this single seed vortex ring manages to start the turbulent 
burst at sufficiently low mutual friction dissipation |22| . 

For the present purpose, to investigate the single-vortex instability, we need an 
injection technique where the seed vortices are initially far apart at low applied cf 
velocity. We can use either remanent vortices [10] or the curved peripheral vortices 
of the equilibrium vortex state, when the sample and rotation axes are not perfectly 
aligned [15] . In both cases rotation is increased rapidly from the initial state at to a 
final constant value Qf where the evolution is recorded at constant external conditions. 
For remanent vortices the initial state is at zero rotation, ft{ = 0, while in the case of 
the equilibrium vortex state fl\ has some constant low value. In a strict sense these two 
methods do not represent injection into flow at constant cf velocity, like the Kelvin- 
Helmholtz instability or the Kibble-Zurek mechanism in neutron irradiation. However, 
in practice they achieve the same result, namely placing evolving vortices in rotating 
cf. In the onset regime, 0.3 T c < T < 0.6 T c , the probability of turbulence depends 
primarily on the final rotation velocity Of and only weakly on the acceleration O used 
to reach Of. We use O ~ 0.02rad/s 2 , which in practice mimics a step increase to J?f. 

Experimental setup: — The measurements are performed in a rotating nuclear 
demagnetization cryostat in which the liquid 3 He sample can be cooled below 0.2 
in rotation up to 3rad/s. The temperature is determined from the frequency shifts in 
the NMR spectra |23l 124] and below 0.3 T c from the damping of a quartz tuning fork 
oscillator [25] . The sample container (Fig.™ is a fuzed quartz tube of radius R = 3 mm 
and length L = 110 mm, filled with liquid ^le at a pressure of P — 29 bar. An aperture 
of 0.75 mm diameter in the center of the bottom end plate restricts the flow of vortices 
into the sample from the heat exchanger volume below. 

The continuous- wave NMR absorption line shape (see insert in Fig. [6| is measured 
non-invasively with two detector coils at both ends of the sample cylinder. The signal is 
recorded with constant frequency excitation by sweeping the polarizing magnetic field. 
The number of vortex lines N in the central vortex cluster is obtained from the NMR 
line shape either experimentally, by comparing to a reference spectrum which has been 
formed with a known number of vortices [22], or theoretically, by comparing to a cal- 
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Fig. 5 (Color online) 3 He sample setup. The oval 2 X 15-turn superconducting NMR tank 
circuit coils are connected inductively via a two-turn pick-up loop to a room temperature 
preamplifier. This provides weak coupling to the pick-up coil, so that the LC resonator is 
not excessively loaded and a Q-value of ~ 6000 is achieved. Two solenoidal superconducting 
magnets |22| (which are not shown) provide the homogeneous axially oriented polarizing fields 
for NMR. 



culated reference, obtained from calculations of the order parameter texture [26] . The 
latter method is applied in this work with an overall accuracy which we estimate to be 
within ±30%. The calibrations are affected by the misalignment between the rotation 
and sample axes. It was measured to be 0.64°. In the following we make use of this 
experimental artifact which breaks cylindrical symmetry and makes the equilibrium 
vortex state at some low initial rotation velocity 0{ ^ useful as a reproducible source 
of seed vortices. 

Such measurements are possible only if the critical velocity v c ,exp(T, P) is stable 
and well-behaved. Since surface defects and dirt on the cylinder wall act as sites for 
nucleation, pinning, and even trapping of vortices, the quartz walls are carefully etched 
and cleaned. In spite of this some variation in critical velocity is observed from one cool 
down to the next, indicating that perhaps frozen gas particles are involved. The sample 
container in Fig. [5] has been in continuous use since a few years, with occasional warm 
ups to liquid nitrogen temperatures to clean the dilution refrigerator circulation from 
air plugs or to room temperature to modify the experiment. During the last 12 months 
~ 90 % of cool downs at 0.8 rad/s to below 0.20 T c remain vortex-free, while at 0.9 rad/s 
only ~ 20% of such attempts are successful. Before that, the same container could be 
regularly cooled down in the vortex-free state at 1.2 rad/s. The exclusion of isolated 
surface defects which trap vortices is no simple task for a large cylinder like that in 
Fig. [5] Note that below 0.3 T c even a single expanding remnant may start a turbulent 
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Fig. 6 (Color online) NMR record of vortex generation in applied cf. (Insert) NMR absorption 
line shapes at 0.51 T c and fl{ = 0.92rad/s. The two line shapes represent: (i) The state after 
increasing rotation from the initial equilibrium vortex state at Q\ = 0.02rad/s to i?f. Here the 
total number of vortices is still close to the initial value TV; ss 5. This state is distinguished 
by the large shifted cf peak on the left, (ii) Final equilibrium vortex state with N eq 840 
(at t > 500 s), which is marked by increased absorption bordering to the Larmor edge on 
the right. Both line shapes have been measured with the top detector. (Main Panel) The cf 
peak height as a function of time after increasing rotation to i?f at 0.50 T c (which is below 
T on = (0.54 ± 0.02) T c ). The reduction in peak height represents the gradual increase in the 
number of vortex lines in the central cluster. The outputs from the top and bottom detectors 
are shown. These are not identical since the vortex formation rates TV and the moments in time 
when the peak heights collapse may differ along the sample. The sawtooth waves represent 
the NMR field sweeps around the location of the cf peak which is shifted far from the Larmor 
value. (Bottom) Rotation drive i?(t) as a function of time. Time t = is placed at the point 
where the final rotation velocity Q{ is reached. 



burst and will then transfer the sample to the equilibrium vortex state. Probably the 
maximum vortex-free flow is here limited by isolated bad spots of dirt or defects on 
the cylindrical wall where vortices can be trapped as small loops, while they try to 
annihilate in zero flow. A trapped loop has a certain critical cf velocity at which it can 
start expanding. This velocity depends on the radius of the loop and its orientation. 
Thus the trapping site controls the flow velocity where the first trapped remnant starts 
evolving. During annihilation the isolated traps are randomly loaded with a remanent 
loop and thus the critical velocity varies from one measurement to the next. 

Measuring procedure: — In the present work the evolving seed vortices are either 
remanent vortices (see Fig. [2| or vortices curving to the cylindrical side wall in the 
equilibrium vortex state (see Fig. |4j. Accordingly, no long temperature sweeps are 
needed here. The rotation drive Q is simply changed at constant temperature according 
to a protocol which for remanent vortices is shown in the inset of Fig. [2] The initial 
state is first formed by decelerating from high rotation with a large number of vortices 
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Fig. 7 (Color online) Number of vortex lines N(t) in central vortex cluster as a function of 
time during continuous vortex formation owing to the single-vortex instability. Two measure- 
ments arc shown which start from equilibrium vortex states at i?j = 0.02 and 0.06rad/s. The 
calculated total number of vortices JV; is initially approximately 5 and 33, respectively, while 
the number of seed vortices A/1, which connect to the cylindrical wall, is 3 — 4 in the former 
and about 20 in the latter case. Owing to difficulties in initial equilibration and the presence 
of a finite annihilation barrier I15| at i7; = 0.02rad/s, the extrapolated number of initial vor- 
tices appears to be 3 - 4 times larger than the calculated estimate at very low velocities like 
0.02 rad/s. The vertical arrows indicate when the turbulent burst occurs in these two examples. 
As in Fig. [6] the burst does not take place inside either of the detector coils, but somewhere 
else in the long column. Since T on increases with A/1, T on has a different value for the two 
cases in this figure: it is 0.538 T c for 17; = 0.02 rad/s, while it is 0.551 T c for Q\ = 0.06 rad/s. 
(Insert) Rotation drive Q(t) and cf peak height measured with the bottom coil at 0.94T on . 
Note that Figs. [3] [6] [7] and |10| - 113| illustrate the single-vortex instability in a long column 
(with length/diameter ~ 20), which is sampled with two detectors, to record the time evo- 
lution simultaneously in two places. Figs. [8] [9] and |14| in turn, describe the situation when 
the column is divided with an A-phase separation layer |17l in two independent samples with 
length/diameter ~ 10, so that the initial equilibration time to the equilibrium vortex state at 
!?i becomes shorter. 



to some low rotation velocity i?;, where rotation is maintained constant for a period 
At. If i?; = 0, then the waiting period At at zero applied flow controls the number of 
remanent vortices [TO], as seen in Fig. [3] If Qi ^ 0, then we generally choose At = 300s, 
which allows the vortex array to approach closer to its equilibrium state. To start 
the single-vortex instability, rotation is next increased by a fixed increment AQ at 
Q ~ 0.02 rad/s 2 to Qf, where it is kept constant and the evolution is recorded. In 
Fig. [6] the NMR response is shown for an example case in the onset regime T w T on , 
where vortex generation starts spontaneously and is finally terminated in a turbulent 
burst. 

The global superfluid cf from the rotation increase to fif produces a large absorption 
peak in the NMR spectrum which is shifted far from the Larmor resonance. In the main 
panel of Fig.[5]the height of this cf peak is monitored. The reduction in peak height as a 
function of time measures the increase in the number of vortices N in the central cluster. 
Well above T on the cf peak height remains constant, as no new vortices are generated, 
but in the onset regime T ~ T on the height may decrease continuously, as seen in 
this example in Fig. [6] The initial slow rate of height reduction we attribute to vortex 
generation by the single- vortex instability. According to the calculated calibrations of 
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Fig. 8 (Color online) Initial rate of vortex generation, TV (at ! = 0), as a function of normal- 
ized temperature T/T on for events with long precursory vortex generation (tj, > 20 s). These 
measurements are performed with the sample divided in independent top and bottom sections, 
separated by a narrow layer of 3 He-A as barrier [2]. Initially both sections are in the equilib- 
rium vortex state at 17;. To start the generation of new vortices, a rapid increase in rotation by 
Afi = 0.7 rad/s is applied. The plot shows that at the final velocity 17f = 17; + AQ events with 
a measurable slow TV are all in the onset regime T fb T on : at high temperatures T/T on > 1.2 no 
cases of dynamic vortex generation were observed, while at low temperatures T/T on < 0.6 all 
measurements ended in a turbulent burst which developed too fast to provide a measurement 
of TV. (Insert) The average TV of the data in the main panel plotted as a function of 17;. Thus 
TV increases approximately linearly with 17;, in other words with the initial number of seed 
vortices A/j . The seed vortices are here the curved peripheral vortices of the equilibrium vortex 
state which at one end connect to the cylindrical side wall. As before, the measured value of 
T on is different for each of the three values of 17; in this figure. Also, vortex generation by the 
single- vortex instability is a stochastic event; this is the origin for the scatter. 



the cf peak heights, the measured TV corresponds to adding a rectilinear vortex line 
every few seconds to the central cluster. The final sudden collapse in height (after about 
50 s in the top and 85 s in the bottom detector) marks the arrival of the equilibrium 
number of vortices to the respective detector coil. 

The collapse of the cf peak is the signal that the turbulent burst has occurred. 
From the site of the burst a vortex front propagates both up and down along the 
rotating column [8] . When the front passes through a detector coil, the cf peak height 
drops to zero. Above 0.4 T c the longitudinal propagation velocity Vp of the front is 
approximately the same as that of the end point of a single vortex while it spirals along 
the cylindrical wall [27], Vhz ~ a(2R (in an originally vortex-free rotating column). 
Recently measurements [28] on the front velocity were extended to temperatures below 
0.2 T c . Using these later values of Vp = a e ffOR and correcting them for the momentary 
number of vortex lines TV in the central cluster around which the front spirals, Vp fa 
a e g v{Q, TV, R), we calculate from the time delay between the collapse of the cf peaks 
in the top and bottom coils the time tb and location z\> of the turbulent burst. In the 
example of Fig. [6] the measured delay of 30 s places the burst at a height zj, = 76 mm 
above the orifice at time % — 32 s (measured from the moment when the rotation drive 
reached Of). 

The analysis of the measured turbulent bursts allows us to conclude that multiple 
bursts, which would occur almost simultaneously, but in different locations along the 



13 




Fig. 9 (Color online) N data from the onset regime in Fig. [s] arranged independently of the 
Qi value in four bins [(0 < TV < Is -1 ), (1,2), (2,3) and, (N > 3s -1 )] in different intervals 
of the normalized temperature T/T on . This histogram of 140 data points characterizes the 
tail of the TV distribution towards slow rates as a function of temperature and illustrates the 
frequency at which such processes occur around T on . Although slow rates (N < 3s — x ) are 
virtually nonexistent below 0.6 T on , this data set does not display a clearly increasing rate 
with decreasing temperature (after averaging over J?i). 



rotating column, are not frequent in the onset regime and clear candidates of such 
events have not been identified. This conclusion is based on the continuous well-behaved 
behavior of the measured data on Vp> and z^. We believe that in the onset regime the 
probability of the turbulent burst is still small and the propagation of the vortex front 
so rapid that it is unlikely for bursts to start at two different locations in sufficiently 
close proximity in time. 

The central characteristic of the single-vortex instability in Fig. [6] is the slowly 
decreasing cf peak height which measures the rate of vortex formation N. A sufficiently 
long period of slow peak height decay for this type of measurement is observed only 
in the onset regime, T ~ T on - At lower temperatures the burst time becomes very 
short after any small rotation increase AQ and our measurement too slow for resolving 
such events. The instability is easier to monitor at lower pressures where longer burst 
times it, are observed [4]. In Fig. [6] T on = 0.54 T c is defined as the average of a series 
of measurements on the transition to turbulence, which fit a normal distribution with 
a half width <tt = 0.02 T c . Here and in all later examples, T on is measured in each 
case separately for the appropriate conditions of that particular measuring situation, 
as was done in Fig. [3] 

To summarize, the collapse of the cf peak height in Fig. [6] is caused by the arrival 
of the vortex front. The front moves with a velocity which depends on the number of 
vortices N in the central cluster at height z, before the front is about to pass at z. 
In an ideal case, where the vortex instability occurs continuously and randomly in the 
sample, one might expect that the site of the burst is randomly distributed along the 
z axis. Varying Q\, we can change in a controlled manner the number of curved seed 
vortices TV; which connect initially to the cylindrical side wall. In addition, by changing 
AQ or temperature, we control respectively the applied flow velocity v(Q, R, N) or the 
damping a(T). By studying the dependence of T on on these variables one can analyze 
how they influence the onset of turbulence (compare Figs.[3]and 14 1. Next we are going 
to focus on the properties of the precursor, represented by the approximately linear 
N(t) in Fig. [£] by examining its characteristics in the onset regime T ~ T on , such as 
the vortex formation rate N, and the distributions of the burst time % and the burst 
locations Zj,. 
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Fig. 10 (Color online) Time t^ needed to start the turbulent burst during continuous vortex 
generation at constant rotation Q{ = Q\ + AQ, plotted as a function of temperature around 
the onset T T on . The data have been collected at five different values of Q\ and three values 
of AQ = 0.7, 0.9, and 1.3 rad/s. Since Ton depends on both Q\ and AQ, in each case the 
appropriate measured value of T on is used for normalizing the temperature axis. The data for 
J?i = come from measurements on remanent vortices with AQ = 0.7rad/s, as in Fig. [3] No 
A-phase barrier field is applied in these measurements and thus the sample is here twice as 
long as in Fig. [8] The plot suggests that on average tj-, increases with normalized temperature 
T/T on and decreases with Q[. 



Experimental results: — Fig.[7]shows two examples, after conversion from cf peak 
height to vortex number N(t). When vortex generation proceeds slowly, the rate N is 
initially of order 1 vortex/s, but increases at later times and becomes more nonlinear, 
until the vortex front passes through the coil and iV jumps to nearly iVeq. As seen 
in Fig. [7[ generally the initial rate N(t = 0) increases with increasing since the 
number of seed vortices A/j, which connect to the cylindrical side wall in the initial 
equilibrium vortex state, increases with 17;. 
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Fig. 11 (Color online) Burst time t^, averaged for the data in Fig. |ip| and analyzed for a 
dependence on the applied flow velocity ~ AQ (left), on the normalized temperature T/T on 
(center), and on the number of seed vortices ~ Qi (right). The averaging has been done by 
dividing the data points in bins as a function of the respective variable and by denoting the 
standard deviation in each bin with uncertainty limits. For instance, on the left we see that t^ 
decreases with increasing flow velocity v(Q,N, R), since this velocity is roughly expressed by 
the rotation increase AQ. Here the solid curve is the average (t^,) = 42/AQ. 
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Fig. 12 (Color online) Temporal distribution of turbulent bursts in the onset temperature 
regime. The burst time data from Fig. |10| is here organized as a histogram for events where 
ib > 20 s is long enough to be resolved. The solid curve is an approximation to the tail of 
the probability distribution for and represents the fitted ex pon ential oc exp (— tb AH/47) j» 
exp (— tb/(it>)j where (ij,) is the average ib a t given Afl (Fig. nT] left). 



More statistics on the precursor are presented in Fig. [8] collected from measure- 
ments similar to those in Fig. [H] The initial rate of vortex generation N(t = 0) is 
compared here in the onset regime T ~ T on for different Q± and thus for different num- 
ber of evolving vortices TV;, keeping the rotation drive ~ Qf — Q t = AQ = 0.7rad/s 
constant. In the inset the average of the measurements is seen to depend roughly lin- 
early on Q\, as was concluded in the context of Fig. [7] The data in Fig. [8] includes 
only events with long burst times > 20 s, where the slope of the cf peak height with 
time can be clearly determined (while events with % < 20 s can be counted, but the 
value of is not resolved) . The data is shown again by the histogram of Fig. 9] but 
now including also the fast events, where t\y < 20 s and N > 3vortices/s. Fig. K) shows 
that, on average, N increases in this data set with decreasing temperature. In addition 
we observe no events with a slow measurable N below 0.6 < T/T on . We thus have to 
conclude that extended precursory vortex generation with a burst time > 20 s can 
only be observed in the onset temperature regime T ~ Ton- 

In Fig. [10] a different data set with measurements on the burst time is examined. 
Since our measurement captures efficiently only events with long burst times, Fig. [10] 
represents the tail (t^ > 20 s) of the burst-time distribution. These 73 data points are 
almost half of all the measured turbulent events in the temperature interval 0.9 < 
T/T ou < 1.1 in the present data set. Thus events with > 20 s and prominent slow 
vortex formation are relatively frequent close to T on . The equations of vortex motion 
roughly scale with the effective rotation drive, which in this case is ~ AQ. As the data 
have been collected with three different values of AQ, we plot the result as t\, AQ 
which corrects for differences in the drive sufficiently well. The main conclusions to be 
drawn from this data are illustrated in the three plots of Fig. |11| where the averages 
are shown as a function of the rotation increase AQ, normalized temperature T/T on , 
and the rotation velocity of the initial equilibrium vortex state J?;. The interpretation 
of these plots is that the burst time ty, decreases rapidly with increasing applied flow 
velocity ~ AQ, increases with temperature T/T on , and decreases with with J?;, i.e. 
the number of seed vortices Mi- In Fig. [l2]the tail of the burst time distribution is 
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Fig. 13 (Color online) Spatial distribution z\ y of turbulent bursts for the data in Fig. |ip[ 
organized as a histogram along the z axis of the long sample in Fig. [5] The bursts appear to 
be randomly distributed, except for a strong preference for the region below the bottom coil. 
This additional mechanism, which enhances the probability of the turbulent burst, we assume 
to be associated with the presence of the orifice below the bottom coil. 



shown as a histogram, indiscriminately for all t\> data in Fig. |10| As seen here, in the 
onset regime, T w T on , the probability for large burst times decreases approximately 
exponentially. 

To summarize we conclude that both data sets, displayed in Figs. [8] and [10] illus- 
trate that well-resolved single-vortex instability events are distinguished by slow vortex 
formation N(t — 0) and a long burst time tb- Such events can be found (i) in the onset 
temperature regime, T w T on , (ii) by starting from a state with a small number of seed 
vortices A/j, and (iii) by keeping the applied flow velocity (~ Ail) as small as possible. 

In Fig. [13] the spatial distribution of the turbulent bursts from Fig. [l0]is displayed 
along the z axis of the sample. As expected, the location z\> of the burst is approximately 
evenly distributed along the column. This supports the notion that the generation of 
new vortices occurs randomly with equal probability along the entire cylinder. The 
exception is a clear preference for the region below the bottom coil. The breakdown 
of these events with zj, < 10 mm according to their J7j values shows that the orifice 
is a large perturbation for small vortex clusters. If J?; = 0.02 rad/s, then the cluster 
radius R ~ R(i7[/ Of)? is approximately equal to the radius of the orifice. Even with 
Q\ — 0.06 rad/s the two might be comparable, since this comparison is affected by 
the centering of the orifice on the bottom plate and the inclination of the cylinder 
and rotation axes. At higher J?; values the cluster apparently covers the orifice more 
efficiently and its perturbing effect fades away. Surprisingly, no cases of turbulent bursts 
are present in Fig. [13] which would have been started by remanent vortices at the orifice 
(i.e. with J?; = 0), but this observation might change with a larger sample of measured 
events. 

Discussion: — We thus find that the presence of the orifice promotes the proba- 
bility of the single vortex instability and concentrates turbulent bursts in the section 
below the bottom detector coil. This is also seen in direct measurements on the indi- 
vidual transition temperatures T on of the top and bottom sections of the long sample 
cylinder. In Fig. [14] the cylinder has been divided in two disconnected halves with a 
narrow transverse layer of 3 He-A, as first discussed in Ref. [17]. The AB interface acts 
as a barrier [2] for vortices and thus the two sections can be sampled separately with 
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Fig. 14 (Color online) Onset temperature of turbulence for the upper and lower sections 
of the sample in Fig. [5] The division into these two sections is described in Ref. |10| . The 
measurements are performed similar to those in Fig. [3] starting from an initial state at zero 
rotation with remanent vortices left over after an annihilation period At = 2 min. Comparing 
the two sections, we recognize that T on is higher in the bottom section, i.e. the probability 
of turbulence at a given temperature is enhanced in the bottom section. The main difference 
between the two sections is the orifice in the lower part. This suggests that the vicinity of the 
orifice is particularly propitious for starting the single-vortex instability. The solid curves are 
normal distributions with a half width cr-p = 0.02 T c , centered around T on . 



the detector coils at each end of the cylinder. Comparing the four T on distributions 
for remanent vortices in Figs. |3]and|14| it is seen that the T on values are all distinctly 
different, but in increasing order form a sequence which can now be explained based 
on the previous conclusions. These four T on distributions have been measured in the 
same set of measurements, without ever warming above dilution refrigerator temper- 
atures. This means that they are closely comparable and representative of the same 
experimental environment. Nevertheless, comparing the two distributions in Fig. |14[ 
it is evident that a relatively large difference of 0.05 T c separates the Ton values of the 
top and bottom sections. This difference is large enough so that the two distributions 
do not overlap. We conclude that the probability of the turbulent burst increases at a 
given temperature when the bottom section is included in the sample — or that the 
vicinity of the orifice is a particularly efficient environment for starting the turbulent 
burst. 

In earlier work 29 it was proposed that a large number of vortices might leak 
simultaneously through the orifice as a sudden burst into the sample volume and that 
this is the origin for the abrupt transition from the vortex-free to the equilibrium 
vortex state. Below the orifice vortices are typically present in the equilibrium vortex 
state in most situations, since two quartz tuning fork oscillators are located there 
25 . These provide ample opportunity for vortex pinning and trapping. Moreover the 
bottom surface in this volume is a rough sintered heat exchanger. However, considering 
Figs.[3]and[l4]together, a massive sudden leakage of vortices through the orifice appears 
now questionable. Comparing the bottom section in Fig. [l4]to the corresponding plot 
for the long sample in Fig. [3] (with At — 2 min on the right), we see that T on is not 
determined by only the orifice, but also depends on the length of the cylinder above 
the orifice: The twice larger sample length L of the long sample causes the onset T on 
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to increase by 0.04 T c . This feature is consistent with the general notion that in the 
long sample there are twice as many seed vortices to begin with, the evolving vortices 
spend more time spiralling along the cylinder wall, and are more likely to suffer the 
instability, than in the short sample. In contrast, an event with lots of vortices leaking 
through the orifice should be insensitive to the length of the sample above the orifice 
and seems less likely to explain the measurements on T on . 

Nevertheless, Fig. [13] shows that a large fraction of the turbulent bursts occur in 
the vicinity of the orifice and all measurements on the bottom section show a higher 
value of Ton by about 0.04 T c than an equivalent measurement on the top section. This 
difference is visible also in Fig. |14| Thus we have to conclude that the single-vortex 
instability is more likely at a given temperature when the orifice is included in the 
sample and that the bottom section therefore has a higher onset temperature than the 
top. It is not clear at this point how the presence of the orifice enhances the probability 
of the single-vortex instability, since it affects both the initial configuration of seed 
vortices and their later dynamics after the rotation increase. Here both geometry as 
well as surface roughness could be important. 

In contrast the lower value of Ton for the top section indicates that an isolated 
cylinder (which in this case is closed off by the AB interface barrier) displays a reduced 
probability for the single-vortex instability to occur, in other words the top section, 
with no obvious defects, is closer to an ideal cylindrical sample. A further measurement 
on only the top section, with a reduced number of initial seed vortices, obtained by 
increasing the annihilation period At from 2min to 20min, reduces T on from 0.44 T c 



(Fig. 14 left) to 0.39 T c . This result is what we expect, based on the examples presented 
above: T on decreases if the number of seed vortices is reduced. 

To summarize, in the onset temperature regime (T « T on ) the single- vortex in- 
stability progresses sufficiently slowly in about one third of the measured events so 
that it can be recorded with our measurement techniques. It functions as the precursor 
mechanism which generates new dynamically evolving vortices, until a localized turbu- 
lent burst between interacting vortices becomes possible in a short section (of length 
~ R) of the rotating column. The instability depends foremost on temperature via 
the dynamic parameter f = (1 — a') /a, which for superfluids is the equivalent of the 
Reynolds number of viscous hydrodynamics, namely the ratio of the inertial and dis- 
sipative forces [2]. Our measurements in the onset regime can be interpreted in terms 
of the probability of a single dynamically evolving vortex to undergo the instability 
and to create a new vortex loop which in turn starts to evolve. At a given temperature 
the probability to achieve bulk turbulence depends on the applied cf velocity at the 
container boundary, on the length of the trajectory over which the vortex end trav- 
els along the boundary, and on the total number of vortices which simultaneously are 
dynamically evolving. 

An earlier explanation of the onset of superfluid turbulence was provided by Klaus 
Schwarz who in 1993 concluded (based on his own work and that of others) that a 
set of several vortex mills is required to start and maintain turbulence in channel flow 
|30J. These vortex mills need to act in parallel and have to be located close to the 
entrance of the flow channel. Our results now show that vortex mills are not necessary 
to start turbulence and that there exists a more fundamental instability mechanism, 
namely the single- vortex instability. In principle, the characterization of this instability 
in Figs. [3] and[6]-[l4]can be compared to simulation calculations, to reconstruct a more 
detailed understanding. A step towards this goal is taken in the next section, where 
the instability mechanism is studied in numerical calculations. 
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3 NUMERICAL CALCULATIONS ON SINGLE- VORTEX 
INSTABILITY 

Numerical method: Our calculations [9| are carried out using the vortex filament 
model introduced by Schwarz |31| . With Biot-Savart integration along all vortex lines 
the superfluid velocity field from vortices is obtained from 



The line integral is taken along the vortices and s(£, t) denotes the location of the 
vortex core at time t, while £ is measured along the arc length of the vortex core. In 
the presence of solid boundaries the total superfluid velocity field, v s = v S)W + Vb, is 
modified by the boundary induced velocity v^. At a plane boundary one can use image 
vortices to satisfy the requirement of zero flow through the boundary, n • v s = 0, where 
n is the unit vector along the surface normal. More generally we obtain vj, = V$ 
by solving the Laplace equation V 2 ^> = combined with the requirement that at the 
boundary n ■ \7<P = —n ■ v s ,u>- No surface pinning or surface friction is included, the 
boundaries are assumed ideal which, as far as known, is not in contradiction with 
measurements. Mutual friction in the bulk superfluid is included using the equation of 
motion for the vortex element at s(£, t) 



where the vector s' = ds/d£ is the local tangent to the vortex at the point s(£,t). For 
the mutual friction parameters a(T,P) and a (T,P) we use the 29 bar data measured 
in Ref. [I3j. 

In practice, the implementation of the Biot-Savart integration is performed adap- 
tively, i.e. the number of discretization points along a vortex in evaluating the Biot- 
Savart integral is increased recurrently, until the required accuracy is obtained. The 
vortex is split along its core into line segments whose length is adjusted such that 
shorter segments are used in places where the vortex is more curved or the counter- 
flow is large (enabling smaller wavelength Kelvin- waves). The smallest segment length 
Z\£ limits the time step At which is used to solve the time development of the tan- 
gle with the classical 4th order Runge-Kutta method. The solution of the Laplace 
equation is obtained by discretizing the potential $ within the cylinder (typical grid 
size, eg. radially, Ar = R/15). The resulting sparse matrix equation is then solved at 
each time step, while the spatial derivatives are approximated with finite differences. 
This means that the continuity equation for the superfluid velocity is not accurately 
satisfied. Nevertheless, this scheme is an improvement over much of the earlier work. 

To solve for Vb, one needs to make sure that vortices meet the boundaries perpen- 
dicularly and that in Eq. |l]) one integrates along vortices which form closed loops, 
as noted by Schwarz |32j . The latter requirement is implemented by extending the 
vortices, which terminate perpendicularly on the boundary, to infinity with straight 
vortex line sections. A vortex reconnection is performed when two vortices approach 
each other closer than the maximum resolution (=A£ ~ R/50 typically, measured 
along the vortex core), provided that the resulting configuration has reduced length 
and represents thus a lower energy state. Generally the maximum resolution has minor 
effect on the results in the intermediate temperature regime. Increased resolution slows 
down the calculations and results in larger numbers of tiny vortex loops which in any 
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Fig. 15 (Color online) (Main ■panel) Calculated build up in cumulative number of vortices 
and their reconnections in a rotating cylinder: (0, solid curve) reconnections in the bulk which 
do not change N, (+1, dashed) surface reconnections which add one new vortex loop, (N) total 
number of vortices, (removed, dash-dotted) small loops which form in reconnections mainly 
close to the cylindrical wall, but which are contracting and are therefore removed, (+1, solid) 
bulk reconnections which add one vortex, (-1, solid) bulk reconnections which remove a vortex, 
(0, dashed) reconnections at the boundary which do not change N. (Insert) Averaged number 
of reconnections per second on the cylindrical boundary and in the bulk which add one vortex. 
The large initial peak in the boundary rate represents the starting burst which is required to 
start vortex formation. Parameters: R = 3mm, L = 10mm, f2 = 0.9rad/s, and T = 0.35T C 
(where a = 0.095 and a' = 0.082 [13]). 



case rapidly disappear owing to the finite mutual friction damping. Nevertheless, a suf- 
ficiently fine resolution is needed to display Kelvin-wave excitations, and towards low 
temperatures the resolution needs to increase rapidly. In solving the Laplace equation 
for the boundary condition, a coarser resolution must be tolerated, to avoid too large 
memory consumption. 

Fig. [2] is an example of how calculations can be used to illustrate and interpret 
measurements. This calculation is performed at 0.38 T c and conserves the number of 
vortices during their evolution, after increasing the rotation velocity from zero to Of. 
Comparing to Figs. [3]and |14| we note that in most experimental cases 0.4 T c is below 
Ton, so that the single- vortex instability would interfere and would lead to turbu- 
lence. The origin of this difference between calculations and experiments has not been 
resolved. To start vortex generation in the calculation, often a specially designed un- 
stable starting configuration is required, which creates a larger number of interacting 
vortices [4]. 

Vortex generation in rotating flow: — In Fig. [15] the generation of vortices is 
studied in a short circular cylinder with diameter comparable to height. Account is kept 
of all reconnection processes as a function of time while the sample is evolving towards 
its final state. Vortex formation is initially started from a single vortex ring which is 
placed in the plane perpendicular to the rotation axis at height 0.2 L slightly off center, 
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50.5 s f = 77.5 s 

Fig. 16 Two snapshots from the calculation in Fig. |l 5) In both cases helically twisted younger 
vortices can be seen on the outer circumference of the cluster and older straighter vortices in 
the center. Primarily surface reconnections contribute to the formation of new vortices at 
t < 100 s which explains the many short loops outside both clusters. 



to break cylindrical symmetry (see Ref. [4]). This is an unstable configuration where 
Kelvin waves of large amplitude immediately form and reconnect at the cylindrical 
wall. The end result is a sudden formation of roughly 30 vortices which have one end 
on the bottom end plate and the other moving in spiral trajectory along the cylindrical 
wall. After this initial burst the spontaneous evolution is followed in the calculation, 
the formation of new vortices is noted, and the reconnections of different type are 
classified. 

In Fig. [16] we see snapshots of vortex configurations after 50 s and 80 s. Recently 
formed vortices are here on the outer circumference in helical configuration, while 
further inside the cluster the vortices are gradually relaxing towards rectilinear lines. 
Outside the cluster closer to the cylindrical wall one can see loops of Kelvin waves, 
small separated loops with both ends of the vortex on the cylindrical wall, and even 
closed vortex rings (lower right corner at t = 50 s). 

Returning to Fig. |15| we note that after the initial burst of the first ~ 30 vortices 
N increases first gradually, but after about 50 s the rate N picks up. During the first 
50 s reconnections in the bulk do not contribute to the generation of new vortices, but 
later such processes also start appearing. Surprisingly however, even during the later 
phase a reconnection of a single vortex at the cylindrical wall, while a Kelvin wave 
expands along this vortex, remains the dominant mechanism of vortex generation. 
This is seen from the fact that the curve for N follows closely that of the successful 
surface reconnections (dashed curve marked as In comparison such surface 

reconnections are few in which a small loop is created, but which later shrinks away, 
for instance, because it is adversely oriented with respect to the azimuthal cf (dash- 
dotted "removed" curve [which actually includes such small loops from both the surface 
and the bulk]). The dashed curve denoted as "0" refers to processes where a closed 
vortex ring from the bulk drifts against the cylindrical wall. Such cases do not change 
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Fig. 17 (Color online) Boundary-mediated vortex formation in pipe flow. The calculation 
shows how a seed vortex travels across the cross section of the pipe after its release. The 
perspective view looks straight into the pipe against the flow, the inner ring is at z = 20 mm 
and the outer at z = 70 mm. Originally at f = the center of the seed vortex was at z = 
closer to the pipe wall on the right. The vortex drifts both along the pipe (towards the viewer) 
and across the flow channel (from right to left). Its center section adopts the curvature of the 
pipe and annihilates (at 45 < t < 50 s). The small loops at both ends survive as independent 
vortices, they are formed to satisfy the boundary condition. Reconnection kinks appear on 
both loops (at t = 50 s), but these do not expand. The two loops reorient themselves with 
respect to the flow and then drift across the flow channel in the opposite direction. In this way 
the number of vortices has increased by one. The repolarization of the two new loops happens 
within the time span 50 - 55 s, as seen in Fig. |18| The present figure shows the start of the 
calculation in Fig. |19| and thus the parameters are here, as well as in Fig. |18| the same as in 
Fig. [19] 

the value of N. They require successful vortex- generating reconnections in the bulk 
and consequently the dashed "0" curve emerges only after the solid curve has 

acquired sufficient slope. In contrast the solid "0" curve represents reconnections in the 
bulk between two different vortices which after the first 40 s rapidly becomes the most 
frequent event. These inter-vortex reconnections do not lead to changes in N and are 
primarily associated with processes occurring between the twisted vortices within the 
bundle, where they help to increase the polarization of the vortices along the rotation 
axis. One might ask whether such bulk reconnections nevertheless emit Kelvin wave 
excitations which then propagate to the boundary and lead to loop formation and 
reconnections there. At present there is no clear evidence of that. 

The insert in Fig. [15] compares the rates of vortex generation from reconnections 
at the wall and in the bulk. The dominant role of wall reconnections is compelling. 
Other similar calculations lead to the same conclusion: The reconnection of a single 
vortex at the cylindrical wall is the most important mechanism for the generation of 
new vortices. This process was illustrated by means of a numerical example in Ref. [1] , 
while the task of Fig.[l5]is to provide quantitative estimates of the relative frequencies 
of successful vortex-generating reconnections at the wall and in the bulk. 

The second important consideration is correspondence with measurement. The ob- 
vious difference between calculation and measurement is the ease with which new 
vortices are generated in experiment below Ton, whereas in Fig. [15] the rate of vortex 



23 



s = 64 




z = 54 



Fig. 18 (Color online) Reorientation of the two new vortex loops in Fig. |17| shown at Is 
intervals. On the bottom (t = 50 s) the two loops have just formed on the back wall of the 
pipe and have the wrong orientation with respect to flow, while at the top (t = 57 s) they have 
changed their orientation and start expanding towards the viewer. The arrows on the loops 
show the orientation of the circulation vector k. The parameters are the same as in Fig. |19| 



generation remains modest. No clearly identifiable turbulent burst can be distinguished 



in Fig. 15 The same calculation at a lower rotation velocity of 0.8 rad/s gives a qualita- 
tively similar result, both with respect to N(t) and the break down in different types of 
reconnections, except that all rates are smaller. After about 150 s both the surface and 
bulk rates turn off simultaneously and vortex generation stops at N « 290 vortices, well 
below the equilibrium number iV C q ~ 570. Clearly in this example no turbulent burst 
takes place, which would boost the vortex number up to JVeq. In Fig. 15 at 0.9 rad/s the 
calculation has been continued to 115 s and N « 400 vortices, where vortex generation 
starts to slow down, again well short of N cq w 650 vortices. The calculations are time 
consuming which limits our possibilities to obtain a more comprehensive understanding 
of their predictions. 

It thus appears as if some mechanism is missing from the calculations in comparison 
to experiment, which makes vortices more unstable and adds to the vortex generation 
rate. The difference is less likely to reside in the bulk than on the cylindrical wall, where 
the condition of an ideal solid boundary should be examined closer. Nevertheless, at low 
vortex density Kelvin-wave formation on a single vortex followed by a reconnection at 
the surface is the only efficient mechanism for generating a new vortex. The problem is 
complex, since Kelvin-wave formation depends on the flow velocity and the orientation 
of the flow with respect to the vortex, which in turn change continuously while the 
vortex moves and the Kelvin wave itself starts to propagate. It appears that with 
respect to the reconnection of Kelvin-wave loops at the wall, the rotating cylinder with 
circular cross section is a particularly stable flow geometry. The next example calculates 
the equivalent of Fig. [15] for linear flow in a circular pipe. Here vortices turn out to be 



24 




Fig. 19 (Color online) Calculated build up in the cumulative number of vortices and recon- 
ncctions in linear pipe flow with circular cross section: (0, solid curve) reconnections in the 
bulk which do not change AT, (+1, dashed) reconnections at the boundary which add one new 
vortex loop, (-1, dash-dotted) very small contracting loops both on the surface and in the bulk 
which approach the resolution limit and are removed, (+1, solid) reconnections in the bulk 
which add one vortex, (-1, solid) reconnections in the bulk between a closed loop and another 
vortex which remove one vortex, (N) total number of vortices, and (0, dashed) reconnections 
at the boundary which do not change N. 



less stable than in rotating flow. The reason is the enhanced role of reconnections at 
the wall. 

Vortex generation in pipe flow: — Technically a measurement with seed-vortex 
injection in vortex-free linear pipe flow of superfluid 3 He-B is a demanding task; so 
far such measurements have not been performed. Nevertheless, we present here calcu- 
lations on a circular straight tube which is initially vortex-free. The calculations are 
performed similar to those above on rotating flow, but by approximating the boundary 
conditions with the faster image vortex techniques. The cf is enforced by imposing on 
the superfluid component flow at constant velocity over the cross section of the pipe. 

Technically such measurements could be set up in the following manner: Suppose 
that both the entrance and the exit of the tube are covered with a superleak which 
prevents the flow of the normal component. The superfluid component is forced into 
motion with a piston acting on a large reservoir in front of the superleak-covered 
entrance of the flow tube. Obviously in a real experiment of this kind large numbers of 
vortices would be created in the superleak. These would continuously flood the tube, 
as long as the flow at constant mass rate is maintained. Such a measurement would 
not be informative about vortex generation and the onset of turbulence. To avoid this 
problem, the tube could be bent to a closed ring, in which the flow is created by 
rotating, in the same way as in the rotating cylinder. 

In our numerical calculation we assume ideal laminar flow through the circular 
straight tube. To start vortex generation, we place one straight vortex line in the flow, 
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t = 300 s t= 370 s 



Fig. 20 Two snapshots from the evolution of the vortex tangle in pipe flow in Fig. |l9| On 
the left the tangle is shown 300 s after the start from one seed vortex and on the right after 
370 s. During their evolution the vortices have drifted in the pipe along the z-axis (dimensions 
are given in mm), where z = was fixed at the middle point of the seed vortex at t = 0. The 
parameters are the same as in Fig. |19| 



which stretches from wall to wall across the tube, slightly tilted from the perpendicular 
plane, to break the symmetry. It turns out that the later evolution of the seed vortex 
is rather insensitive to its initial configuration and that our results do not depend 
materially on how the seed vortex was originally placed in the tube. Experimentally 
seed vortex injection can be achieved by creating vortex rings in applied flow by means 
of the neutron capture reaction in neutron irradiation 22 . 

In Fig. [T7] we examine the trajectory of the seed vortex along and across the flow 
channel. A flat velocity distribution v s = lmm/s is imposed on the superfluid com- 
ponent over the cross section (while the motion of the normal component is clamped 
by the superleaks or by a large kinematic viscosity, v n = 0). The vortex drifts with 
roughly the velocity v s along the pipe downstream (since a <C 1), while it also moves 
transverse across the tube, driven by the dissipative mutual friction force oc av s . The 
consecutive configurations of the vortex are shown in Fig. [17] at 5 s intervals. Owing to 
the boundary condition on the wall, the vortex bows out in the center and mimics the 
curvature of the circular pipe wall, while it traverses across the entire cross section. 
Ultimately its center section, which is now aligned along the pipe wall, annihilates. 
Only one small loop from both ends of the original vortex remains, the vestiges from 
the boundary condition. In this example both of these end loops manage to reorient 
themselves with respect to the flow direction (Fig. 18 I and then start an expanding 
motion in the opposite direction across the flow. A successful reorientation is not al- 
ways the case; often one end loop may simply contract and annihilate. However, in 
Fig. [18] the number of vortices starts to grow continuously from one single seed vortex. 
Here the transverse flight time across the flow is approximately 60s (~ 2R/(av s ). Thus 
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Fig. 21 (Color online) Phase diagram of flow states and dependence of the onset temperature 
of turbulence on flow velocity in linear pipe flow. The calculations are started with a single 
initial vortex, as in Fig. |19| At low velocities below the almost vertical critical boundary 
sustained dynamic vortex generation cannot be maintained at any temperature down to below 
0.3 T c and a state of turbulent flow is not created. At higher velocities above this critical 
boundary sustained turbulence exists at and below the onset temperature plotted on the 
vertical scale. The parameters are the same as in Fig. |19| 



after the first 60s we have two vortices, four vortices after ~ 120s, and after ~ 180s 
N starts increasing more rapidly, as seen in Fig. |19| 

Comparing Figs. |15| and 1 19| it is evident that for pipe flow - in contrast to rotating 
flow - there is no difficulty in starting turbulence from a single seed vortex in the 
simulation calculation. The reason is the difference in flow geometry: As seen in Fig. |17[ 
for pipe flow the generation of new vortices does not depend only on the successful 
expansion of Kelvin-waves, but is aided by the boundary condition which gives rise to 
characteristic end loops. Nevertheless, in both calculations it is the interaction with 
the wall which is responsible for the early phase of vortex generation, in Fig. [19] up 
to about 300 s. After about 360 s the generation and annihilation of loops in surface 
reconnections compensate each other and from here onwards the generation in bulk 
becomes mainly responsible for the production of new vortices. 

Interestingly, also in Fig. [19] reconnections between two different vortices in the 
bulk, which do not directly lead to new vortices, soon dominate over all other processes. 
These bulk reconnections abound, when a sufficiently dense tangle has formed and 
vortices traverse across the flow in both directions, as seen in the two snapshots in 
Fig. |20| This leads to rapidly changing configurations in the tangle. Furthermore, since 
all vortices travel downstream, a turbulent plug is formed which hardly if at all spreads 
upstream (pinning is excluded from our numerical model) . Ultimately, when the vortex 
plug reaches the superleak at the exit of the flow tube, the vortices are annihilated and 
the original state of vortex-free flow reappears. 

It is instructive to note that turbulent plugs have been observed in the flow of 
superfluid 4 He along a straight capillary tube with circular cross section [33]. In these 
measurements temperature fluctuations are registered along the tube in steady state 
flow conditions in some velocity range. The fluctuations are interpreted to arise from 
the flow of turbulent plugs along the tube. The plugs consist of a tangle of quantized 
vortices, they extend over a limited length of the tube, and display relatively sharp 
fronts with the laminar flow sections. However, presumably in the case of 4 He-II the 
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Fig. 22 (Color online) Dependence of the onset temperature of turbulence on the number of 
seed vortices in linear pipe flow. The calculations are started with Afi seed vortices which are 
placed transverse across the pipe at equidistant separation from each other, with their center 
poi nts i n the same cross section of the tube at z = 0. Other conditions are the same as in 
Fig. |19| except for v s = 2 mm/s. At larger values of Afi (above 15 seed vortices) the dependence 
saturates at about T on /T c < 0.55. 



transition from laminar to turbulent flow is not governed by the supply of vortices, 
as discussed here for pipe flow of 3 He-B, but the formation of the plugs is regulated 
by the associated fluctuations in flow velocities, similar to the formation and decay of 
turbulent plugs in viscous pipe flow [53] (where the plugs are formed from classical 
eddies). 

In Figs. [21~| and |22| the calculations from Fig. [19] are carried further, to check for a 
dependence of the onset temperature of turbulence on the flow velocity v s and on the 
number of seed vortices A/;. In Fig.[2l]an almost linear dependence is found for T on , as 
a function of the velocity v s which is imposed on the superfluid component. This line 
represents the onset temperature of the single-vortex instability: At and below this 
critical temperature a state of sustained turbulence is observed to develop. At very 
low velocities (v s < 0.4 mm/s) a region is found where no vortex instability and no 
dynamic generation of new vortices occurs at any temperature down to below 0.3 T c . 
This almost vertical boundary we associate with the critical velocity of turbulence at 
vanishingly small mutual friction. In Fig. [22] in turn, we see that the onset temperature 
depends on the initial number of seed vortices with which the calculation is started: A 



roughly linear dependence T oa (Afi) is found. Figs. 21 and 22 are thus both consistent 
with what we expect on the basis of the rotating experiments (compare to Figs. [8] 
and 10 1. In our calculations on rotating flow these dependences have not emerged 



with equal clarity, owing to the higher stability of evolving vortices in the rotating 
cylinder. In contrast pipe-flow calculations, performed with the same conventions and 
approximations, straightforwardly lead to the expected relations. 



4 SUMMARY 

A quantized vortex is a topologically stable structure of the superfluid order parameter 
field. In principle, it should be possible to account for the appearance of every new vor- 
tex. In 4 He-II this has notoriously been a difficult task: Vortices appear to emerge out 
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of nowhere, without apparent systematics. Various mechanisms have been proposed to 
explain their origin. In 3 He-B at temperatures below 2.3 mK vortex formation processes 
are instabilities and not thermally activated. Vortex formation is here in better control, 
partly owing to the two or three orders of magnitude larger vortex core diameter and 
the reduced influence from surface roughness. Moreover, vortex formation can here be 
examined as a function of a mutual friction dissipation which has a strong, almost 
exponential temperature dependence around the critical regime £ ~ 1. 

With decreasing mutual friction the stability of vortices is reduced and at T ~ 0.6 T c 
turbulence in the bulk volume becomes possible. This hydrodynamic transition from 
laminar to turbulent dynamics takes place at a somewhat higher temperature than 
where an isolated evolving vortex might become unstable in the rotating column. Before 
turbulence can be started from single vortices evolving at low density, new expanding 
loops have to be generated by a precursor mechanism. This happens via the single- 
vortex instability for which the probability rapidly increases with reducing friction. 
Thus the cascade process, the single-vortex instability followed by a localized turbulent 
burst in the bulk volume, becomes possible. Experimentally this is observed as an 
abrupt change in the stability of the dynamics, manifested as a sudden transition to 
turbulence within a narrow temperature interval. The simplest means to investigate 
the instability is to measure at constant rotation velocity the onset temperature T on 
of turbulence after the introduction of a controlled number of seed vortices. Examples 
of such measurements are shown in Figs. [3] and |14[ where T on has been determined 
in each of the four graphs. In Sec. [2] we have analyzed individual data points in such 
graphs, in order to characterize the properties of the single- vortex instability. It thereby 
turns out that in the onset temperature regime, T « T on , a fraction of the transitions 
to turbulence display prolonged precursory vortex formation at slow rate, before the 
turbulent burst in the bulk sets in. At temperatures further below the onset regime 
the instabilities proceed too rapidly to be captured with our measuring techniques. 

Our numerical simulation calculations in Sec. [3] show that interactions of the evolv- 
ing seed vortex with the cylindrical wall in the presence of rotating counterflow is 
the predominant source for new vortices in the low density regime, before interactions 
between vortices in the bulk become possible. In these calculations the walls are repre- 
sented with the boundary conditions of an ideal solid surface. The experimental results 
suggest, however, that surface properties or geometrical features do influence the onset 
temperature, as seen in Fig. |13| More realistic boundary conditions might therefore 
be needed and might reduce the main disagreement between the rotating experiments 
and present calculations, namely enhance the probability of wall reconnections and 
the formation of new expanding loops in the calculations. In contrast, our calculations 
on flow in a straight circular pipe suggest that in this geometry wall interactions lead 
to the generation of new vortices at higher temperatures and at lower flow velocities 
than in rotation. Of these two types of flow, rotation, which promotes polarization of 
vortices along the rotation axis, is a more stable environment for dynamically evolving 
vortices. This is in particular the case for a long smooth-walled circular column with 
good alignment along the rotation axis. The aligned circular column appears to be a 
most stable special case [8] compared to, for instance, a column which is inclined by 
a large amount with respect to the rotation axis or one with square cross section [35]. 
Similar differences in stability are known to apply for viscous pipe flow with circular 
versus square cross section of the flow channel [33] . 

Finally we note that the single-vortex instability starts in the rotating column a 
sequence of events, which have been described in a recent review [8]. Of central interest 
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here is the propagation of polarized vortices along the rotating circular column and how 
this is changed by the increasing turbulent influence with decreasing temperature .28: . 
These phenomena elucidate superfluid turbulence in the T — > limit, which recently 
has been studied also in two other types of measurements, by monitoring the decay of 
homogeneous isotropic turbulence, after the external pumping is switched off [351 136] . 
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